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摘要 : 给 定 连通 地 图 M = (V,E,F)， 在 每 个 节点 v 处 赋 于 一 个 自 旋 ， 每 条 边 上 的 两 个 节点 
相互 作用 ， 得 到 地 图 M 上 的 伊 辛 模型 . 由 伊 辛 模型 的 对 侦 变 换 ， 给 出 了 自 对 偶 地 图 上 伊 
辛 模型 相 变 临界 点 K" 的 上 下 界 : 0.2406 < K" <0.7218 ， 并 且 得 到 三 维 自 对 偶 伊 辛 模型 的 
相 变 临界 点 : K* = 0.5269 . 
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Ising model and critical point of phase transition on self dual map 
ZHANG Xiao-Wu 


Abstract: Given connected map M=(V,E,F), each node V is assigned a spin, two nodes on each 
edge interact, the Ising model on map m is obtained.Dual transformation from Ising model, The 
upper and lower bounds of the critical point K* of the phase transition of Ising model on the self 
dual map are given:0.2406 < K* 0.7218, the critical point of phase transition of the 
three-dimensional self dual Ising model is obtained:K*=0.5269. 
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1 引 


设 有 连通 地 图 ”M =(V,E,)， 对 每 一 个 节点 veV 处 赋 于 自 旋 o,, 每 条 边 


贡 


e= (由 的 两 个 节点 之 间 有 相互 作用 量 /，,，， 如 此 得 到 地 图 M 上 的 伊 六 模型 

当地 图 M =(V,E,F) 是 一 个 环形 图 时 ， 伊 辛 模型 “是 一 维 的 ，M 是 连通 度 
大 于 3 的 平面 图 时 ， 伊 辛 模型 是 二 维 的 ，M 的 亏 格 不 等 于 1 时 ， 伊 辛 模型 是 三 
维 以 上 的 

伊 辛 模型 是 统计 力学 的 一 个 多 体 相互 作用 模型 , 可 以 用 来 描述 非常 广泛 的 物 
理 现象 .如 晶体 的 磁性 , 合金 中 的 有 序 无 序 转变 和 相 变 临界 现象 等 . 一 维 、 二 维 仇 
辛 模型 都 有 精确 解 ， 目 前 为 止 , 三 维 伊 辛 模型 还 没有 得 到 精确 解 "用 低温 展开 、 
高 温 展开 、 重 整 化 群 和 Monte Carlo 模拟 等 方法 ， 能 够 得 到 三 维 伊 辛 模型 的 一 些 
近似 解 & 

参考 文献 中 [5] 介 绍 了 三 维 伊 六 模型 的 数学 结构 与 精确 解 猜 想 ， 从 拓扑 、 代 
数 和 几何 角度 对 三 维 伊 辛 模型 的 数学 结构 进行 了 评述 . 分 析 了 三 维 伊 辛 模型 的 转 


移 矩 阵 、 拓 扑 理论 中 的 纽 结 变换 、Yang-Baxter 方程 和 四 面体 方程 之 间 的 关系 . 
文献 [9] 中 给 出 了 自 对 偶 地 图 和 无 手 征 性 纽 结 的 构造 方法 ， 在 文献 中 [10]， 
也 给 出 了 可 道 纽 结 的 构造 方法 . 纽 结 的 无 手 征 性 就 是 对 应 于 伊 辛 模型 上 的 自 对 
侦 变 换 性 ， 纽 结 的 可 逆 性 就 是 对 应 于 伊 辛 模型 格 点 之 间 相 互 作用 的 对 称 性 . 
文中 利用 地 图 M 上 伊 注 模型 的 对 偶 变 换 ， 给 出 了 自 对 偶 地 图 上 伊 注 模型 相 
变 临 界 点 的 上 下 界 ， 并 且 得 到 三 维 自 对 偶 伊 辛 模 型 的 精确 相 变 临界 点 . 


2 地 图 上 伊 辛 模型 的 配 分 函数 展开 
设 有 连通 地 图 M =(V,E,F)， 其 上 伊 辛 模型 的 Hamiltonian 量 ” 是 : 


H{o.,}=— py Plowo,o, -4pBYo, 


(u,v)eE veV 


其 中 8 是 外 加 磁场 ,地 图 的 每 个 节点 v 是 一 个 蝇 格 , 自 旋 o,=+1, 每 条 边 e = (u,v) 


上 两 个 节点 之 间 的 相互 作用 量 是 ,4=24spys(s+1) 是 磁 窍 ,是 磁 量 子 


数 , = , 这 里 如 是 玻 尔 效 曼 和 常数 ， 了 是 品格 系统 的 温度 . 


设 每 条 边 上 两 个 节点 之 间 的 相互 作用 量 都 是 J/ ， 伊 辛 模型 的 配 分 函数 是 : 
Ov(B,7T)=> >,exp[-H{o,}] =>> [exp[6yoc,] [explBuBo,] 


veV ov Oo, oO, (u,v)eE veV 


其 中 exp[x]=e”. 


， 则 Hamiltonian 量 是 


假设 外 加 磁场 强度 B=0， 天 = 


H{o,}=-K > ao，， 


(uv)eE ' 
由 于 (oo0,) =1 ，exp[Ko,o,]=e’™””"=chK+o,0o,shK =chK(l+o,o,thK) ， 
配 分 函数 为 : 
QV|,T)= > >,exp[-H{o,}= > > [| [ehKQ+o,o,thK)] 
V 


veV ov Oo, oO, (uv)eE 
=(chK) DY TI dto,ow) (w=1hKk). 
oO, Oo, (uv)eE 


将 配 分 函数 的 乘积 项 展开 得 : 
QWV|,T)=(chR) YD tw oo+tw 》 oo,00+.). 


oi=tl cl=tl (i,))eE (Gi,))eE (kK,DeE 


设 n(r) 是 在 给 定 唱 格 上 使 用 了 个 键 所 能 画 出 的 图 形 数 ， 且 这 些 图 形 要 求 每 个 
节点 都 是 由 偶数 个 键 连接 ， 则 配 分 函数 的 高 温 展 开 是 : 


2dy|,7)= ACh) DY nCr)w , (n(0)=D). 


在 基态 下 系统 所 包含 的 所 有 自 旋 都 指向 同一 方向 ， 对 应 总 能 量 = 一 7|E|. 


反 转 一 个 节点 的 自 旋 , 会 减少 4, 个 同 向 最 近邻 对 , 同时 产生 了 4d, 反 向 最 近邻 对 ， 


系统 的 总 能 量 增加 了 24&,7 ,这 里 d, 是 节点 v 的 度数 . 设 |V| 表示 两 个 自 旋 上 -上 


最 近邻 对 的 总 数目 ，|V| 表示 两 个 自 施 下 -下 最 近邻 对 的 总 数目 ， 


V| 表示 两 个 


自 旋 上 -下 最 近邻 对 的 总 数目 ， 并 且 |V| +|V|_+|V|， =| 下 ， 系 统 的 Hamiltonia 
量 为 


HV|,) 一 -7(V|,, 十 VI -|v|, ) 时 -7(E| -2 人 |，) 2 


配 分 函数 的 低温 展开 式 是 : Cdy|,7)= esay mne2*,(m(0) = 


其 中 m(r) 表 示 这 样 的 数目 :设置 蝇 格 中 |V| 个 自 旋 方向 ,， 使 之 有 7 个 反 向 最 近 
邻 对 ，m(7) 表 示 的 就 是 满足 如 此 条 件 之 设置 的 数目 , r=0、d,、… . 
3 伊 辛 模型 配 分 函数 的 对 侦 变 换 


设 有 连通 地 M=(V,E,F) ， 对 偶 地 hM =( 已 天 ) ,EH 


只 = 本 = 由. 


2-25，,M 可 定向 
2- 上 1M 不 可 定向 “ 


M 上 的 晶 格 是 M 上 的 对 偶 唱 格 ， 


-jarlel- zo | 


其 中 加 是 M 的 欧 拉 示 性 数 ，6， 为 M 的 亏 格 数 . 

从 给 定 的 唱 格 M 出 发 ， 在 这 个 唱 格 上 有 一 个 键 数 为 7 的 闭 圈 图 形 ，r 个 键 对 
应 有 (7) 个 闭 圈 , 构建 它 的 对 偶 晶 格 M ， 然 后 在 闭 圈 内 部 的 节点 上 给 定 某 一 方 
向 的 自 旋 ， 外 部 节点 给 定 反 方向 的 自 旋 . 这 样 押 得 的 图 形 是 对 偶 品 格 上 具有 7 个 
反方 向 最 近邻 对 的 财 圈 图 形 ，7 个 反方 向 键 对 应 有 mm，(m) 个 财 圈 . 相反 地 ， 蝇 格 


一 


M 上 茶 图 形 具 有 7 个 反方 向 最 近邻 对 ， 用 路 径 长 度 为 r 的 闭 圈 图 形 来 表示 ，7 个 


反方 向 最 近邻 对 对 应 有 mm (7) 个 闭 圈 图 形 ， 构 建 它 的 对 偶 品 格 M ,这 个 闭 圈 图 


形 束 是 对 偶 晶 格 上 的 键 数 为 7 个 的 闭 图 图 形 , + 个 键 对 应 有 nn.(7) 个 闭 图 , 由 地 图 
的 对 侦 性 得 : ny (7)=m.(7), my (7)=n.(7). 


设 K* = ， 以 使 得 thK* = e2x ， 即 : sh(2K)sh(2K*)=1. 


Qu (WV|,T) =e WY my re =e* Yn.) 六 = . 
r=0 


r=0 


Ci (有 7) 到 2 chk) Dn ) | 


加 
所 以 得 到 配 分 函数 的 对 偶 变换 为 : Ov (V|,T)=2"" (shK'chK*) ?Q,.(|F|,T”). 
4 自 对 偶 伊 辛 模型 的 相 变 临界 点 
在 M =M ， 则 M 是 自 对 偶 的 地 图 ， 其 上 的 唱 格 模型 是 自 对 偶 的 伊 注 模型 ， 


这 时 有 : 
队 |=| 列 = 隐隐 = 加 + 隐 = x ,lBl=2NV|- zy. 


加 
在 临界 状态 时 ，T = 六 ，Ov (| mm=O (yl7)，2 (KecpK ) ? =1. 
_ 固 加 av 中 a 


[el | 1-2— M4. 
2M(shK’chK’) 2 =1, (sh(2K') 7 =2 2 ,sh(2K')=2 =2®, 


| V 
| | | ， 则 相 变 临界 点 下 -mnOvw + V2 TD 


El le 


2 MM 
DE 


对 于 连通 地 图 MM ， 


vi-1g|als 3WIVI- 


2|y| V| 4 Te | | 
1<1 | <1， val eR li) 
局 可 E Pi 所 有 n+ . n(2+V5) 


自 对 偶 地 图 M 上 伊 辛 模型 的 相 变 临界 点 KK 满足 : 0.2406 < K* <0.7218. 
对 于 三 维 伊 辛 模型 ， 节 点 数 |V|=mnl ， 度数 d,=6 ， 边 数 |BI=3mnl, 设 M 


在 曲面 上 的 谍 入 是 自 对 偶 的 地 图 ， 面 数 |[F|= mnl, xw =-1121 ，M 是 六 正则 地 


图 M 可 定向 时 , 亏 格 数 名 ,= 了 (2+ wn) =1+ =， 其 中 mn 为 偶数. M 不 可 定向 
时 ，é& =2+nml. 


于 是 有 : K， 一 有 -mnGE+WA4+D=05269， 
根据 现 有 文献 [5-8]， 由 Monte Carlo 模拟 和 重 整 化 群 等 近似 方法 得 到 相 变 
临界 点 的 近似 值 是 : 


=0.2174,0.22200,0.2216544,0.2217,0.3074 


六 


ksTc 


三 维 自 对 侦 伊 辛 模型 有 较 高 的 对 称 性 和 稳定 性 ， 其 相 变 临界 点 KK" = 0.5269 比 模 
拟 得 到 的 临界 点 都 要 大 ， 文 献 [5] 中 由 二 个 猜想 得 的 的 Curie 温度 相 变 临界 点 ， 
V5-1 


和 地 存在 于 黄金 点 exp( 2K.)=— =0.61800339887, 所 得 到 的 


SE 
uy 
er 


K. = 0.2406059 与 文中 的 临界 点 下 界 精确 值 5mG 十 = 0.2406059 近似 相等 ， 


两 个 临界 点 的 误差 很 小 . 
5 结论 

根据 地 图 上 伊 辛 模型 的 对 偶 变 换 和 地 图 是 自 对 偶 的 充 要 条 件 , 本 文 得 到 了 自 
对 偶 地 图 的 上 伊 辛 模型 相 变 临界 点 ”的 上 下 界 : 0.2406< 开 "<0.7218 ， 并 且 给 
出 三 维 自 对 偶 伊 辛 模 型 的 相 变 临界 点 : K* = 0.5269 . 
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